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ABSTRACT 

On 2011 July 14, a new magnetar candidate, Swift J1822.3— 1606, was identified via a rate trigger 
on the Swift /Burst Alert Telescope. Here we present an initial analysis of the X-ray properties of the 
source, using data from the Rossi X-ray Timing Explorer, Swift, and the Chandra X-ray Observatory, 
spanning 2011 July 16 - September 22. We measure a precise spin period of P = 8.43771963(5) s and a 
spin-down rate of P = 2.97(28) x 10~ 13 , at MJD 55761.0, corresponding to an inferred surface dipole 
magnetic field strength of B = 5.1 x 10 13 G, the second lowest thus far measured for a magnetar, 
though similar to IE 2259+586 as well as to several high-magnetic field radio pulsars. We show that 
the pulsed X-ray flux decay in the 2-10 keV band is best fit by an exponential with a time constant of 
16.4±0.3days. After increasing from ~35% during the first week after the onset of the outburst, the 
pulsed fraction in the 2 - lOkeV band remained constant at ~45%. We argue that these properties 
confirm this source to be a new member of the class of objects known as magnetars. 
Subject headings: stars: magnetars — pulsars: individual (Swift J1822. 3-1606) — stars: neutron — 
X-rays: bursts — X-rays: general 



1. INTRODUCTION 

On 2011 July 14, the Swift/Burst Alert Telescope 
(BAT) triggered on several bursts of hard X-ray emission 
from the direction of a previously unk nown source, sub- 
seque ntly named Swift J1822.3— 1606 (jCummings et all 
1 2 1 lh . The spin period of the new source wa s soon iden- 
tified to be P = 8.4377s (jGogus et al.ll201lD and obser- 
vations with Swift X-ray Telescope (XRT) were used to 
furth er localize the position of the source (|Pagani et al.1 
l20Tl . 

The simplest interpretation of these initial ob- 
servations is that Swift J1822.3— 1606 is a magne- 
tar. However, because of suggested similarities be- 
tween Swift J1822.3— 1606 and the unusual Be X- 
ray binary J1626— 5156, namely, a similar pulsed frac- 
tion and IR counterpart, it was suggested that per- 
haps Swift J1822 . 3— 16 6 is also a Be X-ray b i- 
nary (Gog us et al.1 120111: iBandvopadhvav et all |2011[) . 



This hypothesis was challenged by Irlalpernl ( 2011 ). 
who, due to the lack of an optica l counterpart 
(|de Ugarte Postigo fe Munoz-Dariasll201lD . argued that 
the properties of Swift J1822.3— 1606 are in line with 
those of other magnetars, in particular, the transient 
magnetar XTE J1810— 197. Here we report on timing 
and flux properties of Swift J1822.3— 1606 that rule out 
a Be binary, and instead agree strongly with the magne- 
tar identification. 

Magnetars are neutron s tars powered by the decay 
of extreme magnetic fields ([Thompson fc Duncan! [19951 : 
iThompson fc Duncanf! 9961; IThompson et al.H2002D . Soft 
Gamma-ray Repeaters (SGRs) and Anomalous X-ray 
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Pulsars (AXPs) are thought to be observational manifes- 
tations of magnetars, and we refer to them as such here. 
Magnetars exhibit X-ray bursts, spin periods (P) in the 
relatively narrow rangeQ of 2 - 12 s, and are observed 
to spin down with time. As inferred from P and spin- 
down rate P via the conventional vacuum dipole brak- 
ing model that implies dipolar magnetic field strength 
B = 3.2 x 10 19 (PF) 1/2 G, their inferred dipole magnetic 
fields tend to be significantly larger than for classical pul- 
sars, typically B > 10 14 G. Magnetars are also observed 
to display variability in nearly all of their observed prop- 
erties: flux and pulsed flux variations, spectral varia- 
tions, pu lse profile variab i lity, timing noise and glitches 
(e.g.. see lMereghetti|[2008t IRea fc Espositoll2011l for re- 
views). 

The true underlying distribution of magnetic field 
strengths of magnetars is an open question. The recent 
discovery of a magnetar with a very low inferre d magnetic 
field of B < 7.5 x 10 12 G (SGR 0418+5729; IRea et al.l 
2010), has raised the question of how low the dipolar 
field can be for magnetar-like activity to be observed. 
Moreover, the discovery of a magnetar-like outburst from 
the high-g rotation-powered pulsar (RPP) J 1846-0258 
(jGavriil et ai][200l Kumar fc Safj-Harbl200l ) highlights 
the question of whether rotation-powered pulsars and 
magnetars represent distinct classes of object, or if mag- 
netars represent the high-i? tail of a single population. 

In this Letter, we present analyses of Rossi X-ray Tim- 
ing Explorer (RXTE), Swift, and Chandra X-ray Ob- 
servatory (CXO) data. We perform a phase-coherent 
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timing analysis and show that the spin evolution is con- 
sistent with a constant spin-down rate. We present an 
analysis of the decay of the source's pulsed intensity af- 
ter the outburst, and show that the pulsed fraction is 
steady after increasing during the first week after the 
outburst. We argue that these properties confirm that 
Swift J1822.3— 1606 is a new Galactic magnetar. 

2. OBSERVATIONS 
2.1. Swift Observations 

On 2011 July 14 at 12:47:47 UT the Swift/BAT trig- 
gered on a rate increase from the previously unknow n 
source, Swift J1822.3-1606 (|Cummings et al.l 120111 ). 
The XRT began observing Swift J1822.3-1606 on 2011 
July 15 at 17:43:20 UT. A total of 24 observations were 
taken with the XRT (target IDs 32033 and 32051). The 
observations span 67 days: 2011 July 15 - September 20 
(MJD 55758 - 55824), for a total exposure time of 46 ks. 

Cleaned Swift data in both windowed-timing (WT) 
mode with 1.8-ms time resolution and photon-counting 
(PC) mode with 2.5-s time resolution were downloaded 
from the Swift quick-look data archive^. For the lone 
PC mode observation, an annular source region with an 
outer radius of 20 pixels and an inner radius of 5 pixels 
was extracted. The inner radius was excluded because of 
pile-up. For the WT mode observations, a 16 pixel strip 
centered on the source was used. Events in the 2-10 keV 
energy range were then extracted from the reduced WT 
mode data for the timing analysis. The PC mode data 
have insufficient time resolution to be useful for timing. 

2.2. RXTE Observations 

Swift J1822.3— 1606 data were obtained with the P ro- 
portional Counter Array (PCA; I Jahoda et al.l 120061 ) on 
board RXTE. The PCA consists of five xenon/methane 
proportional counter units (PCUs) sensitive to 2 - 60keV 
photons, with an effective area of ~ 6500 cm 2 and a 
- 1° FWHM field of view. We downloaded 23 pub- 
lic RXTE observations of Swift J1822.3-1606 from the 
HEASARC archiveO (observing program P96048). Data 
were collected in GoodXenon mode, which records the ar- 
rival time (1-fJS resolution) and energy (256-channel res- 
olution) of each event. The observations span 62 days: 
2011 July 16 - September 16 (MJD 55758 - 55820), for 
a total exposure time of 119ks. 

For our timing analysis, we extracted 2-10 keV pho- 
tons (because this produced high significance pulse pro- 
files for individuals observations) from the top xenon 
layer of each PCU to reduce contamination from the par- 
ticle background. Data from individual PCUs were then 
merged. When more than one observation occurred in 
a 24-hr period, wc combined the data to produce better 
profiles. 

2.3. Chandra Observations 

CXO ACIS-S observations were made on 2011 July 27 
(MJD 55769.2, ObsID 12612), August 4 (MJD 55777.1, 
ObsID 12613), and September 18 (MJD 55822.7, ObsID 
12614) with net exposures of 15.0 ks, 13.7 ks, and 10.0 ks, 
respectively. All three observations were obtained in 

3 http:/ /swift. gsfc.nasa.gov/cgi-bin/sdc/ql 
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continuous clocking mode, with 2.85-ms time resolution. 
We carried out the data reduction using CIAO 4.3 with 
CALDB 4.4.33. We reprocessed the data to retain events 
on the chip node boundary, then extracted the source 
counts from a 6"-width aperture and the 0.5-8 keV en- 
ergy range for our timing analysis, and 2-8 keV for the 
pulsed fraction analysis. 

3. ANALYSIS & RESULTS 
3.1. Phase- coherent Timing Analysis 

Our timing analysis of Swift J1822.3— 1606 follows 
the common phase-coherent approach, in which we ac- 
cou nt for each rotation of th e pulsar (see, for exam- 
ple iManchesteL^SiHs^IIIZS) an d utilized a combina- 
tion of timing data from RXTE, Swift, and Chandra. 
For data from each telescope, events were reduced to 
bary centric dynamical time (TDB) at the solar system 
barycenter using the XRT p osition of RA=18 h 22 m 18 s , 
Dec=-16°04'26.8" (J2000; iPagani et al1l20H and the 
JPL DE200 solar system ephemeris. Events were then 
binned into time series with 31.25-ms time resolution. 

Each RXTE (2 - lOkeV), Swift (2 - lOkeV), and CXO 
(0.5 - 8keV) time series was folded with 64 phase bins 
using a frequency determined from a periodogram anal- 
ysis. After finding an initial phase-coherent timing solu- 
tion, we used this ephemeris to re-fold all the profiles to 
produce higher quality pulse Times Of Arrival (TOAs). 
For RXTE data, we found that using 128 phase bins 
created good pulse profiles with optimal TOA uncertain- 
ties. For Swift and CXO data, 64-bin profiles resulted in 
higher significance pulse profiles and good TOA uncer- 
tainties. 128- and 64-bin template profiles were created 
by aligning and summing all RXTE profiles. 

To account for the different energy ranges and tele- 
scope responses between RXTE, Swift, and CXO, which 
result in small but significant differences in the profiles, 
we fit for a constant phase offset between TOAs obtained 
from different telescopes. We used a subset of the avail- 
able data to fit for the phase offsets (MJD 55758 - 55781), 
which were then held fixed for the remainder of the tim- 
ing analysis. 

We searched for time variability in the pulse profile 
and found that individual profiles are consistent with the 
template in each case except for the first observation after 
the outburst (96048-02-01-00). In this case, the variation 
is subtle, at the ~0.015 phase level. We accounted for 
this by multiplying the corresponding TOA uncertainty 
by three. 

For each profile, we implemented a Fourier domain fil- 
ter by using six harmonics in the cross-correlation pro- 
cedure with the appropriate 64- or 128-bin template. 
Cross-correlation produces a TOA for each observation 
with a typical uncertainty of ~27ms (0.32% of P) for 
RXTE TOAs, ~59ms (0.70% of P) for Swift TOAs, and 
~23ms (0.27% of P) for CXO TOAs. TOAs were fitted 
with the software package TEMPIlQ; parameters from this 
fit are given in Table [1] 

The top panel of Figure [1] shows timing residuals with 
only spin-frequency (y) fitted, while the bottom panel 
shows residuals with the frequency derivative (y) also 

5 http://cxc.harvard.edu/ciao4. 3 
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fitted. Fitting v improves the RMS residuals by a factor 
of 1.8 and improves the x 2 ' v jv from 3.66/44 to 1.14/43. 
Resulting phase residuals show minimal evidence for an 
unmodelled trend and fitting further parameters does not 
significantly improve the fit. The fitted values of v and v 
imply a surface dipolar magnetic field of B = 5.1 x 10 13 G. 

3.2. Pulsed Count-rate Evolution 

To measure the source intensity and decay after the 
outburst, we first checked for the presence of dust scatter- 
ing rings, such as t hose observed around IE 1547—5408 
(jTiengo et al.ll20"To() . because these will bias the inferred 
source count rate. We checked for deviations from the 
point-spread functions (PSFs) in the first Swift PC mode 
and CXO observations following the outburst. The ra- 
dial profiles were well fit by the model PSFs, hence we 
found no evidence for dust scattering rings. 

To measure the pulsed count rate for each Swift and 
CXO observation, the barycentered time series were 
folded with the ephemeris (see Table Q] and Section 13. ip , 
with 16 phase bins. The pulsed count rate for each folded 
profile in the 2 - lOkeV energy range was then deter- 
mined using a root-mean-squared (RMS) method. We 
calculate the RMS puls ed count rate us ing Fourier com- 
ponents as described in lDib et al.l (|2008D . using 5 Fourier 
harmonics of the pulse profile. The pulsed fraction was 
then determined by dividing the pulsed count rate by the 
total source count rate. 

To measure the pulsed count rate for RXTE data, we 
created phase-resolved spectra with 16 bins using the 
software program f asebirQto fold barycentered photons 
with the ephemeris (Table [T]). We selected photons in the 
2-10 keV energy range and from the top xenon detection 
layer. To minimize instrumental effects, we included only 
data from PCU2. PCUO and PCU1 have lost propane 
layers, adversely affecting their particle background lev- 
elso, and data available from PCU3 and PCU4 were min- 
imal for this source. The RMS method with 5 harmonics 
was then used to measure the pulsed count rate. Because 
the PCA is not a focusing instrument the pulsed fraction 
could not be accurately determined from these data. 

Figure [2] shows the evolution of the pulsed source in- 
tensity (top panel, as measured with Swift and RXTE), 
and the pulsed fraction (bottom panel, as measured with 
Swift and CXO). We present here the pulsed source in- 
tensity evolution only in terms of count rate, instead 
of fluxes. A detailed spectral analysis will follow in a 
subsequent paper. To describe this pulsed intensity de- 
cay quantitatively, we fit the pulsed count-rate evolution 
of Swift J1822.3-1606 in the 2 - lOkeV energy range 
with an exponential decay. Pulsed count rates from a 
given source for RXTE and Swift arc different owing 
to instrumental effects (e.g. effective area), but should 
be offset only by a constant, and trends should be the 
same between instruments (as is observed here). Thus, 
RXTE PCA pulsed count rates were scaled to the Swift 
XRT pulsed count rates. The scaling factor was cho- 
sen to be that which minimized the residuals of the fit 
to the pulsed count rates. The exponential decay is de- 
scribed by F(t) = F p exp-(*- to )/ r +F q where F p is the 
peak count rate, F q is the count rate in quiescence, to is 

7 http:/ /heasarc. nasa.gov/docs/xte/recipes/fasebin. html 
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the time of the BAT trigger in MJD and r is the decay 
timescale in days. With a xt °f 2.69 for 41 degrees of 
freedom, the best-fit decay timescale is 16.4 ± 0.3 days. 
A power-law model was also fit to the decay, but it pro- 
vided a much worse fit with a xt of 79 for 41 degrees of 
freedom. 

The bottom panel of Figure [2] shows the evolution of 
the 2-10 keV pulsed fraction of Swift J1822.3-1606 
for Swift and CXO data. The pulsed fraction appears 
to have increased from ~35% immediately after the out- 
burst until about one week after the BAT trigger. It then 
remained constant at ^45%. 

4. DISCUSSION 

The regular spin-down we report for 
Swift J1822.3— 1606, together with a pulsed flux 
decay similar to that seen in other magnetar outbursts, 
demonstrates that this source indeed shares critical 
properties with other known magnetars; we hence 
classify it as such. The inferred surface dipole B for 
Swift J1822.3-1606 of 5.1 x 10 13 G is smaller than those 
of all but one confirmed magnetar, SCR 0418+5729 
(|Rea et al.ll2010f) . though it is ve ry close to that of AXP 
IE 2259+586 (B = 5.9 x 10 13 G: iGavriil fc Kaspill200i_h . 
As for the latter, Swift J1822.3-1606's B is similar 
to that measured for several RPPs, including the lone 
identified magnetically active R PP, PSR J1846-025 8 
which has B = 4.9 x 10 13 G (|Gotthelf et al.l [20001. 
Although enhanced spin-down could be present if the 
source suffered a , large glitch at the outburst (e.g. 
iKaspi et al.1 120031 ). the duration of our observations of 
Swift J1822.3— 1606 (64 days) significantly exceeds the 
glitch recovery time scales measured in previous large 
magnetar glitches, so any contamination is likely to be 
small. 

Interestingly, the most recently discovered magnetars 
(including SGR 0418+5729 and Swift J1822. 3-1606, 
but also SGRs 0501+4516, 1833-0832 and 
Swift J1834.9-0846; Gogiis, et al. 2010a,b; Rea et al. 
2010, Kuiper & Hermsen 2011) all have B < 1 x 10 14 G, 
effectively lowering the 'average' inferred B-field for 
magnetars. This raises the question of what is the true 
distribution of magnetar field strengths. In Figure [31 we 
plot the distributions of the periods and inferred surface 
dipolar magnetic fields of all confirmed magnetars, as 
well as of all known B > 10 12 G radio pulsars. 

The P distribution of magnetars remains narrow, span- 
ning less than an order of magnitude, in great contrast to 
those of radio pulsars. We note that the period distribu- 
tion for magnetars is statistically consistent with being 
flat; a fit to the histogram mean yields a \ 2 V °f 0.63. 
This is somewhat surprising and needs to be addressed 
in any future population studies. The B distribution for 
observed magnetars, by contrast, appears more peaked, 
even with the relatively low B we have measured for 
Swift J1822.3— 1606, albeit with the possible low tail rep- 
resented by SGR 0418+5729. If the magnetars represent 
the high-S tail of the B distribution of all non-recycled 
neutron stars, then their peaked B distribution may indi- 
cate that burst rates of these sources fall rapidly with de- 
creasing B, since the primary mechanism for their discov- 
ery is through their bursting behavior, particularly with 
the Swift BAT and Fermi GBM. This wo uld be broadly 
consistent with the theoretical findings of iPerna fc Pons! 
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(2011i), who calculate expected burst rates in neutron 
stars and show that all neutron stars can show magnetar- 
like bursts, though the burst rate drops with decreasing 
B and increasing age. That the distribution of radio 
pulsars plummets rapidly above B ~ 10 13 G could gen- 
uinely be due to fall-off in the intrinsic distribution, or 
to smaller radio beams for these preferentially longer- 
period pulsars, but may also indicate that radio emission 
is har der to produce at higher B (e.g. iBaring fe Harding! 
120011 ) . Although radio emission has been observed for 
three mag netars (jCamilo et all I2006L 120071 : E evin et aLl 
I2010T ) . their radio spectra and radiative stabilities are 
vastly different than for conventional radio pulsars, sug- 
gesting a distinct emission mechanism. Hence, if indeed 
we can 'unify' radio pulsars and magnetars as described 
above, we expect the lowest B magnetars to be the like- 
liest to produce observable traditional radio-pulsar-like 
radio emission, and the highest B RPPs to be the like- 
liest to show magnetar-like X-ray outbursts. The lat- 
ter is consistent wi th the behavior of PSR J1846— 0258 
(|Gavriil et al.l l2008). while as yet, there is no evidence to 
support the former prediction. Such a unification would 
also suggest that there exists a large population of as yet 
undetected B ~ f 13 G neutron stars that could be one 
day detected via their occasional magnetar-like bursts. 

The characteristic age (t c = P/2P) for 
Swift J1822.3-1606 is 450 kyr, comparable to that 
of AXP IE 2259+586 but much larger than for all 
other magnetars, except for SGR 0418+5729 (>24Myr). 
IE 2259+586 is firmly associated with the supernova 
remnant C TB 109 which is thought to have an age 
of ^9 kyr ([Sasaki et al.l 120041 ). much smaller than the 
AXP's characteristic age of 250 kyr. This demonstrates 
that t c for magnetars can be unreliable true age in- 
dicators, which perhaps can be ascribed to epochs of 
greater spin-down torque early in the star's history, 
when the dipolar field was larger. Indeed this hypoth- 
esis can also explain the large t c of SGR 0418+5729, 
although the latter's greater distance from the Galactic 
Plane suggests it may indeed be the oldest known 
magnetar ( but see lAlpar et alj 120111 for an alternative 
discussion). iTurolla et al.1 ( 20111 ) suggest that the 
observed properties of SGR 0418+5729 are consistent 
with it being an aged magnetar in which the external 
B has decayed significantly. Further evidence for 
SGR 0418+5729 being older than other magnetars 
includes lower energy bursts and a very low quiescent 
luminosity. Whether Swift J1822.3— 1606 more closely 
resembles SGR 0418+5729 or IE 2259+586 remains to 
be seen. The measurement of its flux in quiescence and 
the discovery of an associated supernova remnant would 
help resolve this issue. 

The post-outburst pulsed count rate evolution of 
Swift J1822.3— 1606 is best characterized by an expo- 
nential decay with timescale of 16.4+0.3 days. This 
is similar to what has been obse rved after sev- 
eral other magnetar outbursts (e.g. iRea et al.1 120091 : 
IGotthelf fc Haloernl [20051 : IGavriil et al. l I2008D. though 
power-law decays are also common (e.g. iWoods et al.l 
120011: iKouveliotou et al.|[200l . iLvubarskv et al.l ([2002D 
predicted power-law decays in magnetar outbursts as- 
suming crustal cooling following an impu lsive heat in- 
jection. iBeloborodov fc Thompson] (|2007l ). by contrast, 



considered magnetar outbursts as a result of sudden 
twisting of magnetic field lines in the magnetosphere, 
with the relaxation a result of their untwisting. In this 
model, the relaxation is predicted to be approximately 
linear, a functio nal form observ ed in one outburst o f 
IE 1048.1-5937 (IDib et all I207)9l) . IBeloborodov! (pOOflh 
showed that other functional forms for the decay, includ- 
ing exponential, are possible, as the untwisting is pre- 
dicted to be strongly non-uniform, being erased by a 
propagating 'front' whose speed depends on the initial 
twist configuration. The diversity in functional forms for 
the flux decays of magnetars hence favors the model put 
forth bv iBeloborodovr ([2009T) . 
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Figure 1. X-ray timing residuals for Swift J1822.3— 1606 spanning 2011 July 16 - September 22. RXTE points are shown as filled red 
circles, Swift points are black crosses, and CXO points are blue triangles. The top panel shows residuals with only v fitted, while the 
bottom panel shows timing residuals with v also fitted. Parameters from this fit arc given in Table [l] 
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Figure 2. Top: Pulsed count rate evolution of Swift J1822.3— 1606 in the 2 - lOkeV energy range following its outburst. Swift data are 
shown as black crosses, RXTE data arc filled red circles. The dashed line indicates the best-fit exponential decay with r = 16.4 ± 0.3 days. 
Bottom: Pulsed fraction measurements in the 2-10 keV with Swift data as black crosses and CXO data as blue triangles. 



Table 1 

Spin Parameters for Swift J1822.3-1606. 



Parameter Value 



Dates (Modified Julian Day) 55758 - 55822 

Epoch (Modified Julian Day) 55761.0 

Number of TOAs - RXTE 21 

Number of TOAs - Swift 22 

Number of TOAs - CXO 3 

v (Hz) 0.1185154336(8) 

v (Hz 2 ) -4.2(4) X 10~ 15 

RMS residuals (ms) 23.4 

RMS residuals (periods) 0.0028 



Derived Parameters 

P (s) 8.43771963(5) 

P 2.97(28) x 10~ ] 

Surface Dipolar Magnetic Field, B (G) 5.1(2) X 10 13 

Spin-down Luminosity, E (ergs -1 ) 2.0(2) X 10 31 

Characteristic age, r c (kyr) 450(40) 



Note. - 



Uncertainties arc formal lcr TEMPO errors. 



The Spin-down of Swift J1822.3-1606 
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Figure 3. (Top) Histograms of the spin periods of observed magnetars (for which P has also been measured; see the online magnetar 
catalog, http://physics.mcgill.ca/~magnetar/main.html; solid) and known radio pulsars having inferred B > 10 12 G (from the ATNF 
pulsar catalog at http://www.atnt.csiro.au/research/pulsar/psrcat/ dotted). Swift J1822.3— 1606 is indicated. Note the different vertical 
axes for the radio pulsars (lett) and magnetars (right). (Bottom) Histograms of the inferred B-field strength of the same magnetars (solid) 
and the same radio pulsars (dotted). Swift J1822.3— 1606 is indicated, as is the upper limit for SGR 0418+5729 (Rca ct al. 2010). Note 
again the different vertical axes. 



